function [ Fgrid, Bgrid ] = Matrix_Bivariate_4point( grid , omega, iter )
%-----------------------------------------------------------------
% Input:
% grid(m,n,:,:) - a grid over Z^2 of matrices
% omega - tension parameter
% iter - # of iteration needed
% Output:
% Fgrid(2^(iter+1)*m-5(1+2+4+...+2^iter),2^(iter+1)*m-5(1+2+4+...+2^iter),:,:)
%-----------------------------------------------------------------
% ABSTRACT
% Tensor product. the middle point is insert in the order y->x of
% calculation
%-----------------------------------------------------------------
% NIr Sharon, 21-05-12
%-----------------------------------------------------------------
Fgrid = grid;
Bgrid = grid;

for i=1:iter
   [Fgrid] = Matrix_Bivariate_4pointRefinement( Fgrid , omega );
   
   [Bgrid] = Matrix_Bivariate_4pointRefinement_boundary_V2( Bgrid , omega );
end


end

